Dephasing in quantum chaotic transport: a semiclassical approach 
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We investigate the effect of dephasing/decoherence on quantum transport through open chaotic 
ballistic conductors in the semiclassical limit of small Fermi wavelength to system size ratio, 
Xf/L <g 1. We use the trajectory-based semiclassical theory to study a two-terminal chaotic dot 
with decoherence originating from: (i) an external closed quantum chaotic environment, (ii) a clas- 
sical source of noise, (iii) a voltage probe, i.e. an additional current-conserving terminal. We focus 
on the pure dephasing regime, where the coupling to the external source of dephasing is so weak 
that it does not induce energy relaxation. In addition to the universal algebraic suppression of weak 
localization, we find an exponential suppression of weak-localization oc exp[— t/t^,], with the dephas- 
ing rate t^ 1 . The parameter f depends strongly on the source of dephasing. For a voltage probe, 
f is of order the Ehrenfest time oc 1ii[L/Af]. In contrast, for a chaotic environment or a classical 
source of noise, it has the correlation length £ of the coupling/ noise potential replacing the Fermi 
wavelength Ap. We explicitly show that the Fano factor for shot noise is unaffected by decoherence. 
We connect these results to earlier works on dephasing due to electron-electron interactions, and 
numerically confirm our findings. 

PACS numbers: 05.45.Mt,74.40.+k,73.23.-b,03.65.Yz 



I. INTRODUCTION 

A. Dephasing in the universal regime 

Electronic systems in the mesoscopic regime are 
ideal testing-grounds for investigating the quantum-to- 
classical transition from a microscopic coherent world 
(where quantum interference effects prevail) to a macro- 
scopic classical world 1 -. On one hand, their size is inter- 
mediate between macroscopic and microscopic (atomic) 
systems, on the other hand, today's experimental control 
over their design and precision of measurement allows one 
to investigate them in regimes ranging from almost fully 
coherent to purely classical 2 ^! 4 -. The extent to which 
quantum coherence is preserved in these systems is usu- 
ally determined by the ratio /t c \ of the dephasing time 
T0 to some relevant classical time scale t c \. For instance, 
t c i can be the traversal time through one arm of a two- 
path interfcrometer^ii, or the average dwell time spent 
inside a quantum dot 8 -' 9 - ' 10 ' 11 ' 12 . In a given experimental 
set-up, can often be tuned from > t c \ (quantum 
coherent regime) to <C t c \ (purely classical regime) by 
varying externally applied voltages or the temperature of 
the sample. 

Coherent effects abound in mesoscopic physics, the 
most important of them being weak-localization, uni- 
versal conductance fluctuations and Aharonov-Bohm 
interferences in transport, as well as persistent cur- 
rentsSi 3 -^. The disappearance of these effects as dephas- 
ing processes are turned on has raised lots of theoret- 
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menta l 26 i 27 ' 28 i 29 ' 30 i 31 ' 32 ' 33 interest. Focusing on transport 



through ballistic systems, dephasing is usually investi- 
gated using mostly phenomcnological models of dephas- 
in g 16 i 17 i 18 i 19 i 20 i 21 , the most successful of which are the 
voltage-probe and dcphasing-lead model a 16 ' 17 . In these 
models, a cavity is connected to two external, L (left) and 
R (right) transport leads, carrying Ar and Al transport 
channels respectively. Dephasing is modeled by connect- 
ing a third "fictitious" lead to the system, with a volt- 
age set such that no current flows through it on aver- 
age. Electrons leaving the system through this third lead 
are thus re-injected at some later time, with a random- 
ized phase (and randomized momentum). These models 
of dephasing present the significant advantage that the 
standard scattering approach to transport can be applied 
as in fully coherent systems, once it is properly extended 
to account for the presence of the third lead. 

Using random matrix theory (RMT), the volt- 
age / dephasing probe models^ 4 , predict an algebraic sup- 
pression of the weak-localization contribution to the con- 
ductance (in units of 2e 2 /h)^, 



wl 

9rmt 



N K N L 



[N R + N L ] 



[1 + T D /% 



(1) 



where — ArAl/[Ar + Al] 2 is the weak-localization cor- 
rection in the absence of dephasing. Similarly, universal 
conductance fluctuations become 1 - 1 - 
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time reversal symmetry. For a fictitious lead connected 
to a two-dimensional cavity (a lateral quantum dot) 
via a point contact of transparency p and carrying A3 
channels one has = mA/hpN^, with the electron mass 
m and the area A of the cavity. Similarly, the dwell time 
through the cavity is given by td = mA/H(Ni, + Nr). 

The dcphasing and voltage probe models account for 
dephasing at the phcnomcnological level only without 
reference to the microscopic processes leading to dephas- 
ing. At sufficiently low temperature, it is accepted that 
the dephasing arises dominantly from electronic interac- 
tions, which, in diffusive systems, can be well modeled by 
a classical noise potential^. Remarkably enough, this 
approach reproduces the RMT results of Eqs. ([TJ) and 
([2]) with set by the noise power. These results are 
moreover quite robust in diffusive systems. They are es- 
sentially insensitive to most noise-spectrum details, and 
hold for various sources of noise such as electron-electron 
and electron-phonon interactions, or external microwave 
fields. For this reason it is often assumed that dephasing 
is system-independent, and exhibits a character of uni- 
versality well described by the RMT of transport applied 
to the dephasing/voltage probe models^. 



B. Departure from RMT universality 

According to the Bohigas-Giannoni-Schmit surmise 3 ^, 
closed chaotic systems exhibit statistical properties of 
hermitian RMT— in the short wavelength limit. Open- 
ing up the system, transport properties derive from the 
corresponding scattering matrix, which is determined by 
both the Hamiltonian of the closed system and its cou- 
pling to external leads^. It has been shown that for not 
too strong coupling, and when the Hamiltonian matrix 
of the closed system belongs to one of the Gaussian en- 
sembles of random hermitian matrices, the corresponding 
scattering matrix is an element of one of the circular en- 
sembles of unitary random matrices 3 ^. One thus expects 
that, in the semiclassical limit of large ratio L/Ap of the 
system size to Fermi wavelength, transport properties of 
quantum chaotic ballistic systems are well described by 
the RMT of transport. This surmise has recently been 
verified semiclassically^. 

The regime of validity of RMT is generally bounded 
by the existence of finite time scales, however, and it 
was noticed by Aleiner and Larkin that, while the de- 
phasing time t$ gives the long time cut-off for quantum 
interferences, a new Ehrenfest time scale appears in quan- 
tum chaotic system in the deep semiclassical limit, which 
determines the short-time onset of these interferences 1 ^. 
The Ehrenfest time te corresponds to the time it takes 
for the underlying chaotic classical dynamics to stretch 
an initially localized wavepacket to a macroscopic, clas- 
sical length scale. In open cavities, the latter can be 
either the system size L or the width W of the open- 
ing to the leads. Accordingly, one can define the closed 
cavity, Tg = A _1 1ii[L/Af], and the open cavity Ehren- 



fest time, t° p = A" 1 ln[W 2 /X F L]^.. The emergence 
of a finite te strongly affects quantum effects in trans- 
port, and recent analytical and numerical investigations 
of quantum chaotic systems have shown that weak local- 
izatio n 14 ' 43 i 44 ' 45 i 46 and shot nois o 47 ' 48 i 49 ' 50 are exponen- 
tially suppressed oc exp[— te/td] in the absence of de- 
phasing (t0 — * 00). Interestingly enough, the deep semi- 
classical limit of finite te sees the emergence of a quan- 
titatively dissimilar behavior of weak-localization and 
quantum parametric conductance fluctuations, the latter 
exhibiting no TE-dependcncc in the absence of dephas- 
ingjjiif&Si^. These results are not captured by RMT, 
instead one has to rely on quasiclassical approache a 14 i 44 
or semiclassical method s 43 ' 45 ! 46 ' 55 ' 56 to derive them. 



C. Dephasing in the deep semiclassical limit 

The behavior of quantum corrections to transport at 
finite te in the presence of dephasing was briefly inves- 
tigated analytically in Ref. [lj], for a model of classi- 
cal noise with large angle scattering, and numerically in 
Ref. [23[ , for the dephasing lead model with a tunnel bar- 
rier. Intriguingly enough, the two approaches delivered 
the same result, that quantum effects are exponentially 
suppressed cx exp[— te/t^]. This suggested that dephas- 
ing retains a character of universality even in the deep 
semiclassical limit. More recent investigations have how- 
ever showed that at finite Ehrenfest time, Eq. ((T|) be- 



comes^ (see also Ref. [24 
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with a strongly system-dependent time scale f . Ref. [25[ 
showed that, for the dephasing lead model, f = t e ' + (1 — 
/o)r E p in terms of the transparency p of the contacts to 
the leads, which provides theoretical understanding for 
the numerical findings of Ref. [23| . If however one con- 
siders a system-environment model, where the environ- 
ment is mimicked by electrons in a nearby closed quan- 
tum chaotic dot, one has f = t£, where 



(4) 



in terms of the correlation length £ of the inter-dot inter- 
action potential. 

On the experimental front, an exponential suppression 
oc exp[— T/T c ] of weak-localization with temperature has 
been reported in Ref. [2{|. Taking ~ T" 1 as for de- 
phasing by electronic interactions in two-dimensional dif- 
fusive systems, this result was interpreted as the first ex- 
perimental confirmation of Eq. ([3]) . There is no other the- 
ory for such an exponential behavior of weak-localization, 
however, the temperature range over which this experi- 
ment has been performed makes it unclear whether the 
ballisti o 13 ' 15 , ~ T~ 2 , or the diffusive dephasing time 
determines the Ehrenfest time dependence of dephasing 
(see the discussion in Ref. [24[). 
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Weak localization 


Conductance fluctuations 


Shot noise 


System with environment 


T = Te 




no dephasing 


Classical noise (microwave, etc) 
e-e interactions within system 
System-gate e-e interactions 


f = Tc 
f = Tc 


f ~ 0, Ref. [24] 
f ~ Ref. [24] 
t ~ T£ , follows from R,ef . [2^| 


Dephasing lead: 

no tunnel barrier 

low transparency barrier 


T — T E 

op | cl r, cl 

T = Tjf +T B ~ 2t e 


f = 

f ~ 2-rJ; 1 (numerics), Ref. [231 


no dephasing, Ref. [57] 



Table I: Summary of the known results to date on the nature of the exponential term exp[— t/t$\ in the dephasing [cf. Eq. ([3])]. 
Results that are not referenced are obtained in the present article and in Ref. [25[. Here we list the value of f for different 
transport quantities and different sources of dephasing, all in the pure dephasing regime (the phase-breaking regime of Ref. [BH]). 
The parameter £ differ slightly from system to system (see text for details) however it is always related to the correlation length 
of the interaction with the environment. The results of Ref. [24| neglect r^-contributions (so "0" could indicate a f of order 
tj), they are also only valid for t e 1 ^> t$. 



D. Outline of this article 

In the present article, we amplify on Ref. [25| and ex- 
tend the analytical derivation of Eq. (J3]) briefly presented 
there. We investigate three different models of dephasing 
and show that the suppression of weak-localization cor- 
rections to the conductance is strongly model-dependent. 
First, we consider an external environment modeled by 
a capacitively coupled, closed quantum dot. We restrict 
ourselves to the regime of pure dephasing, where the en- 
vironment does not alter the classical dynamics of the 
system. Second, we discuss dephasing by a classical noise 
field. Third, following Ref. [56|, we provide a semiclassi- 
cal treatment of transport in the dephasing lead model. 
For these three models, we reproduce Eq. and derive 
the exact dependence of f on microscopic details of the 
models considered. All our results are summarized in 
Table |H 

The outline of this article goes as follows. In Section 
im we present the treatment of the system-environment 
model, focusing in particular on the construction of a 
new scattering approach to transport that incorporates 
the coupling to external degrees of freedom. We ap- 
ply this formalism to a model of an open quantum 
dot coupled to a second, closed quantum dot. We 
present a detailed calculation of the Drude conductance 
and the weak-localization correction, including coherent- 
backscattering, which explicitly preserves the unitarity 
of the S'-matrix, and hence current conservation. This 
calculation is completed by a derivation of the Fano fac- 
tor, showing that, in the pure dephasing limit, shot noise 
is insensitive to dephasing to leading order. In Sec- 
tion [Till we present a model of dephasing via a classical 
noise field (such as microwave noise). We consider clas- 
sical Johnson-Nyquist noise models of dephasing due to 
electron-electron interactions within the system, and de- 
phasing due to charge fluctuations on nearby gates. In 
Section IIVI we present a trajectory based semiclassical 
calculation of conductance in the dephasing lead model, 



both for fully transparent barriers and tunnel barriers. 
We also comment on dephasing in multiprobe configura- 
tions. Finally, Section [V] is devoted to numerical simu- 
lations confirming our analytical results. Summary and 
conclusions are presented in Section IVI1 while technical 
details are presented in the Appendix. 



II. TRANSPORT THEORY FOR A SYSTEM 
WITH ENVIRONMENT 

In the scattering approach to transport, the system is 
assumed fully coherent and all dissipative processes oc- 
cur in the leads^. Apart from its coupling to the leads, 
the system is isolated. Here we extend this formalism 
to include coupling to external degrees of freedom in the 
spirit of the standard theory of dccohcrcncc. The cou- 
pling to an environment can induce dephasing and re- 
laxation. Here, we restrict ourselves to pure dephasing, 
where the system-environment coupling does not induce 
energy, nor momentum relaxation in the system. In semi- 
classical language, we assume that classical trajectories 
supporting the electron dynamics are not modified by 
this coupling. 

The starting point of the standard theory of decoher- 
ence is the total density matrix r] ioi that includes both 
system and environmental degrees of freedom^. The ob- 
served properties of the system alone are contained in 
the reduced density matrix ?7sys, obtained from rjtot by 
tracing over the environmental degrees of freedom. This 
procedure is probability conserving, Tr [?7sys] = 1, but 
it renders the time-evolution of 77 sys non-unitary, and in 
particular, the off-diagonal elements of rj sya decay with 
time. This can be quantified by the basis independent 
purity^, < Tr [?7s ys ] < 1, which remains equal to one 
only in the absence of environment. We generalize this 
standard approach to the transport problem. 
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Figure 1: Schematic of the system-environment model. The 
system is an open quantum dot that is coupled to an environ- 
ment in the shape of a second, closed dot. 



A. The scattering formalism in the presence of an 
environment 

We consider two capacitively coupled chaotic cavities 
as sketched in Fig.Q] The first one is the system (sys), an 
open, two-dimensional quantum dot, ideally connected to 
two external leads. The second one is a closed quantum 
dot, which plays the role of an environment (env). The 
two dots are capacitively coupled, and in particular, they 
do not exchange particles. Thus current through the sys- 
tem is conserved. We require that the size of the contacts 
between the open system and the leads is much smaller 
than the perimeter of the system cavity but is still semi- 
classically large, so that the number of transport chan- 
nels satisfies 1 <C iVL.R *C L/Xp. This ensures that the 
chaotic dynamics inside the dot has enough time to de- 
velop, Atd ^> 1, with the classical Lyapunov exponent 
A. Electrons in the leads do not interact with the sec- 
ond dot. Few-electron double-dot systems have recently 
been the focus of intense experimental efforts^. Parallel 
geometries, of direct relevance to the present work, have 
been investigated in Refs. [62|l63l |. 

The total system is described by the following Hamil- 
tonian 



be tuned by applying external backgate voltages on a 
given system. 

In the standard scattering approach, the transport 
properties of the system derive from its (TVl, + N-r) x 
(TVl + -ZVr) scattering matrh*i£ 



S = 



SLL s R l 
SLR Sr,r 



(G) 



which we write in terms of transmission (t = slr) and 
reflection (r = s QjQ ,, a E {L, R}) matrices. From S, 
the system's dimensionlcss conductance (conductance in 
units of 2e 2 /h) is given by 



g = Tr(ttt). 



(7) 



To include coupling to an environment in the scattering 
approach, we need to define an extended scattering ma- 
trix § that includes the external degrees of freedom. This 
is formally done in Appendix [5] and our starting point is 
Eq. (|A11[) for the case of an initial product density matrix 

»7tot = Vsys ® ?7env, with Tjiyl = \n)(n\, n € {1, • • ■ , N L }, 
and ?7 G nv, the initial density matrix of the environment. 
We define the conductance matrix 



(m Tr cnv U 
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'/sys 
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(8) 



where Tr onv stands for the trace over the environmental 
degrees of freedom. From this matrix, the dimensionlcss 
conductance is then given by, 



* = EE 

9nm 3 

n— m— 



(9) 



and Eq. (jAllj) reads 



E 

n£R-.m£L 



dqdqodqj, (q |?7env|qo) 



xS r , 



1 (q,q )(s„ m (q, q„) 



(10) 



Eq. fjlOjl is the generalization of the Laudauer-Buttikcr 
formula in the presence of an external environment. It 
constitutes the backbone of our trajectory-based semi- 
classical theory of dephasing. 



It. 



(5) 



Inside each cavity the chaotic dynamics is generated by 
the corresponding onc-particlc Hamiltonian H sys cnv . We 
only specify that the capacitive coupling potential XL is 
a smooth function of the distance between the particles. 
It is characterized by its magnitude U and its correlation 
length £ such that its typical gradient is U/£. Physi- 
cally, £ is determined by the electrostatic environment of 
the system, such as electric charges on the gates defining 
the dots and the amount of depletion of the electrostatic 
confinement potential between the gates and the inver- 
sion layer in semiconductor heterostructures. Generally 
speaking, U and £ are independent parameters and can 
have different values in different systems, and might even 



B. Drude conductance 

The scmiclassical derivation of the one particle scatter- 
ing matrix has become standar d 64 ' 65 i 66 . Once we intro- 
duce the environment we deal with a bipartite problem, 
here we use the two-particle semiclassical propagator de- 
veloped in the framework of entanglement and decoher- 
The extended scattering matrix elements can 
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,67,68 



be written as, 

Smn(q,qo) : 



, (m\y) (y \n) 
dt 1 dV ° l R dV (ar*)H~-*>/» 
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where we take a d sys -dimcnsional. one-particle system 
(throughout what follows we take d sys = 2), and a d- 
dimensional, iV-particlc environment. At this point, § 
depends on the coordinates of the environment and is 
given by a sum over pairs of classical trajectories, labeled 
7 for the system and T for the environment. The classical 
paths 7 and T connect yo (on a cross-section of lead L) 
and qo (anywhere in the volume occupied by the environ- 
ment) to y (on a cross-section of lead R) and q (anywhere 
in the volume occupied by the environment) in the time 
t = ^ = tr- For an environment of TV particles in d 
dimensions, q is a Nd component vector. In the regime 
of pure dephasing, these paths are solely determined by 
H sys and ff e nv- Each pair of paths gives a contribution 
weighted by the square root A^At of the inverse deter- 
minant of the stability matri x 69 ' 70 , and oscillating with 
one-particle (S" 7 and «Sr, which include Maslov indices) 

and two-particle (§ 7 ,r = Jo^^b^ 7 ")' c lr( T )]) action in- 
tegrals accumulated along 7 and T. 

We insert Eq. (fTTj) in Eq. (fTU)) . sum over chan- 
nel indices with the semiclassical approximatio n 45 ' 56 
Yin 1 ' (Vo\ri) (n\y' ) ~ S(y' - y ). For th e environme nt 
we make the random matrix ansatz that (qo|?7cnv|qd) ~ 
(2nH) Nd Z-i v 5(q' - q ), where 5 onv is the environ- 
ment phase-space volume. The dimensionless conduc- 
tance then reads 



9 = 



(2-Kh)- 1 



dtdt' / dq dq / dy / dy 

"onv JO ienv JL JR 

7,r;7'T 



A T At- e *(*v+*-»+*u). (12) 



This is a quadruple sum over classical paths of the system 
(7 and 7'; going from yo to y) and the environment (r 
and T', going from qo to q) with action phases, 

$ sys = [&y (y,y ;t) - Sy (y,y ;t')]/h, (13a) 
$cnv = [S r (q,qo;*)-Sr'(q,qo;*')]M> (13b) 

$u = [S 7 ,r(z/,yo;q,qo;*) - §y,r'(j/,j/o;q ) qo;* / )]/' i - 

(13c) 

We are interested in quantities averaged over variations 
in the energy or cavity shapes. For most sets of paths the 
phase of a given contribution will oscillate wildly with 
these variations, so the contribution averages to zero. In 
the semiclassical limit, Eq. (|12l) is thus dominated by 
terms which satisfy a stationary phase condition (SPC), 
i.e. where the variation of <& sys + $ env + $u has to be 
minimized. In the regime of pure dephasing, individual 
variations of 'I'sys, $cnv and $u are uncorrelated. They 
are moreover dominated by variations of <I> sys and $ env , 
on which we therefore enforce two independent SPC's. 

The dominant contributions that survive averaging are 
the diagonal ones. They give the Drude conductance. 
Indeed setting 7 = 7' and r = V straightforwardly satis- 
fies SPC's over $ sys and "^env These two SPC's require 
t = t' and lead to an exact cancellation of all the phases 



$ sys = $ cnv = $u = 0. The dimensionless Drude con- 
ductance is given by 



at 



dq dq / dy a / dy 

'L JR 



A 2 A 2 
^ (2nh) S CE 



(14) 



From here on, the calculation proceeds along the lines 
of Ref . [1^ . The main idea is to relate semiclassical am- 
plitudes with classical probabilities. This is done by the 
introduction of two sum rules that express the ergodic 
properties of open cavities, Eq. (|15a[) . and of closed ones, 
Eq. (HE}, 
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d0 o d0P 8ys (Y,Y o ;i) [•••] Yo , (15a) 



£>r["']r = / dp dp P env (Q,Q ; *)[■•■ ] Qo . (15b) 
r ^ 

Here, P sys (Y ,Y ;t) = p F cos6 x P sy8 (Y; Y ;t), and 
P sys (Y,Y ;i) and P e nv(Q,Qoi*) are the classical prob- 
ability densities. For the system, we need to take into 
account the fact that particles are injected, which is why 
the classical probability density must be multiplied with 
the initial system momentum p-p cos 9q along the injection 
lead 6 ^. The phase points Y = (yo,#o) and Y = (y,6) 
are at the boundary between the system and the leads. 
In contrast Qo = (qo,Po) and Q = (q, p) are inside 
the closed environment cavity. The momenta are inte- 
grated over the entire environment phase-space, while 
Pcnv(Q 7 Qo! t) will always contain a <5-function which re- 
stricts the final energy to equal the initial one, (i.e. |p| = 
I Po| if all N environment particles have the same mass). 

The average of P sys over an ensemble of systems or over 
energy gives a smooth function. For a chaotic system we 
write 



Psys(Y; Y ; t) 



cos 9 



2 (W L + Wr) t d 



,-t/TD 



(16a) 



Likewise, the average of P env gives 



Penv(Q;Qo;0) = ' (lpl : l , Pol) ; d6b) 



dPol 
^env 



where sjjnv is the size of the hypersurface in the envi- 
ronment's phase-space defined by |p| = |po| (for d = 

2, sinv = (27rpp nv ri cn v) Ar where fl cnv is the area of 
the environment). Inserting Eqs. (|15a[) . (|15b[) . (|16ap 
and (fTobl) into Eq. ([T4]) . we can perform all integrals us- 
ing Jcnv d Qo = Jc„v d qo d PO = Senv and J cnv dQ5(|p| - 

|po|) = sin"' ■ Then, since 7Vl,r = feWi^R/i": we recover 
the classical Drude conductance, 



9 Nt + N, 



R 



(17) 
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C. Overview of the effect of environment on 
weak-localization 

The leading-order weak-localization corrections to the 
conductance were identified in Refs. [IlUzil as those 
arising from trajectories that are paired almost every- 
where except in the vicinity of an encounter. An ex- 
ample of such a trajectory is shown in Fig. [21 At the 
encounter, one of the trajectories (7') intersects itself, 
while the other one (7) avoids the crossing. Thus, they 
travel along the loop they form in opposite directions. 
For chaotic ballistic systems in the semiclassical limit, 
only pairs of trajectories with small crossing angle e con- 
tribute significantly to weak-localization. In this case, 
trajectories remain correlated for some time on both sides 
of the encounter, the correlated region indicated in pink 
in Fig. [2] In other words, the smallness of e requires 
two minimal times, TL(e) to form a loop, and T\y(e) in 
order for the legs to separate before escaping into dif- 
ferent leads. In the case of an hyperbolic dynamics one 
estimates^, 



r L (e) 

Tw(e) 



A _1 ln[e- 
-'He-^W/L) 2 ] 



,-2l 



A 



(18a) 
(18b) 



As long as the system-environment coupling does not 
generate energy and/or momentum relaxation, the pres- 
ence of an environment does not significantly change this 
picture. However it does lead to dephasing via the accu- 
mulation of uncorrelated action phases, mostly along the 
loop, when 7 and 7' are more than a distance £ apart. 
This is illustrated in Fig. [2] for the case when £ is less 
than W. We define a new timescale as twice the time 
between the encounter and the start of the dephasing, 



T^nX-'He-WLf 



(19) 



Dephasing occurs mostly in the loop part. However if 
£ < eL and < 0, dephasing starts before the paths 
reach the encounter. We discuss this point in more detail 
below in Section Hi Gl 



D. Calculating the effect of the environment on 
weak-localization 



In the absence of dephasing each weak-localization 
contribution accumulates a phase difference <5$ S ys = 
E-pe 2 / (A/t) 55 i 71 . In the presence of an environment, an 
additional action phase difference 5&u is accumulated. 
Incorporating this additional phase into the calculation 
of weak-localization does not require significant depar- 
ture from the theory at U = 0. We extend the theory of 
Ref. [HI to account for this additional phase. 

We follow the same route as for the Drude conduc- 
tance, but now consider the pairs of paths described in 
Section III CI above, and shown in Fig. [2l while the en- 
vironment is still treated within the diagonal approxi- 
mation, T' = T. The sum rule of Eq. (|15b|) still ap- 
plies. Though the sum over system paths is restricted 




Chaotic system 
(environment not shown) 

Figure 2: (color online) A semiclassical contribution to weak- 
localization for the system-environment model. The paths are 
paired everywhere except at the encounter, where one path 
crosses itself at angle e, while the other one does not (going 
the opposite way around the loop). Here we show £ > eL, 
so the dephasing (dotted path segment) starts in the loop 
(T 4 > 0). 



to paths with an encounter, we can still write this sum 
in the form given in Eq. (|15a[) . provided the probability 
P sy s(Y,Yo;i) is restricted to paths which cross them- 
selves. To ensure this we write 

P sys (Y,Y ;i) =p F cos6» / dR 2 dRiP sys (Y, R 2 ; t - t 2 ) 



x P sys (R 2 ,R 1 ;t 2 -ti)P sys (R 1 ,Y ;ii).(20) 

Here, we use R = (r, </>), </> € [— 7r,7r] for phase-space 
points inside the cavity, while Y lies on the lead as be- 
fore. We then restrict the probabilities inside the in- 
tegral to trajectories which cross themselves at phase- 
space positions Ri,2 with the first (second) visit to the 
crossing occurring at time t\ (i 2 ). We can write dR 2 = 
Vp sined£idt 2 de and set R 2 = (i"i,0i ± e). Then, the 
weak-localization correction to the dimcnsionlcss conduc- 
tance in the presence of an environment is given by, 



.9"' = 
with, 



JdY J deRe[e^](F(Y ,e)),(21a) 



F(Y ,e) = 2v F sme 



dt 

Tl+Tw ^T l + - 



dh 



t 2 -T L 



dtj 



x PF cos 6» / dY / dRiP sy8 (Y,R 2 ;t-t 2 ) 
Jr Jc 

x Psys(R 2 , Rl! i 2 — il)Psys(Rl, Yo; t\) 

^S^Penv(Q,Qo;*)exp[W* u ]. (21b) 
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Comparison with Eq. (34) of Ref . [45[ shows that the ef- 
fect of the environment is entirely contained in the last 
line of Eq. (|21bj) . At the level of the diagonal approxi- 
mation for the environment, T' = T, one has 

5$ U = \J dr[U(r 7 (r),qr(r)) - ll(iy (r), q r (r))] , 

(21c) 

where r 7 (r) and qr^) parametrize the trajectories of 
the system and of the environment respectively. We note 
that in the absence of coupling, 5&u = 0, the integral 
over the environment is one, and we recover the weak- 
localization correction of an isolated system (cf. Eq. (35) 
of Ref. [H). 

To evaluate Eqs. (|2"Tj) . we need the average effect of the 
environment on the system, after one has traced out the 
environment. For a single measurement, this average is 
an integral over all classical paths followed by the envi- 
ronment, starting from its initial state at the beginning 
of the measurement. We therefore also average over an 
ensemble of initial environment states, or an ensemble 
of environment Hamiltonians, which corresponds to per- 
forming many measurements. For compactness we define 
(■ ■ • )env as this integral over environment paths and the 
ensemble averaging over the environment, 



dQdQ 



(Penv(Q;Qo;*)[---] Qo >. (22) 



Without loss of generality we assume that for all r, the 
interaction ll(r, q) is zero upon averaging over all q. We 
can ensure an arbitrary interaction fulfills this condition 
by moving any constant term in IX into the system Hamil- 
tonian (these terms do not lead to dephasing). Since the 
environment is ergodic, we have 



<lt(r 7 (<),qr(i))) c 



= 0. 



(23) 



Now we use the chaotic nature of the environment 
to give the properties of the correlation function 

(U(r T (t),qrW)lt(ry(t'),qr(t'))> env - 

(i) Correlation functions typically decay exponentially 
fast with time in chaotic systems, with a typical 
decay time related to the Lyapunov exponent^. 
The precise functional form J(A on v|i' — t\) of the 
temporal decay of the coupling correlator depends 
on details of H cnv and IX, however for all practi- 
cal purposes, it is sufficient to know that it de- 
cays fast, and we approximate it by a <5-function, 
J(\ em \t'-t\)~\-i v 6(t). 

(ii) We argue that the spatial correlations of II for two 
different system paths at the same t also decay, 
because the averaging over many paths, and many 
initial environment states act like an average over q. 
We define K ( | r ' — r | /£) as the functional form of this 
decay of spatial correlations, with K(0) = 1. The 
precise form of K(x) depends on details of -ff sys , 



iJ env and 11. In particular, the typical length £ 
of this decay is of order the scale on which ll(r, q) 
changes between its maximum and minimum value. 

Given these arguments, we have 

(U(r 7 (0,qrWMMO,qr(O)) cnv 
= ffitf(|r r (t)- 1^)1/0 5(t-t>). (24) 



Then, 
(<5$ 2 



U/cnv 



= Wo dT2dTl 



x ( [IX (r 7 (r 2 ) , q r (r 2 ) ) - IX (ry (r 2 ) , q r (r 2 ) ) ] 
x [IX (r 7 (n ) , q r (n ) ) - IX (ry (n ) , q r (n )) ] 



= 2^ /V[l-Jf(|r y (T)-r 7 (T)|/0] . (25) 

^env Jo 

We further make the following step-function approxi- 
mation for K, 



K{x) = 9(1 - x), 



(26) 



where the Eulcr 9-function is one (zero) for positive (neg- 
ative) arguments. In principle, this is unjustifiable for 
x ~ 1, however since the paths diverge exponentially 
from each other, the time during which x ~ 1 is of order 
A -1 , while dephasing happens on a timescale which 
is typically of order the dwell time, tq. Thus the step- 
function approximation of K(x) will have corrections of 
order (Atd) -1 -C 1, which we therefore neglect. Once 
we have made the approximation in Eq. (|26[) . we see that 
non-zero contributions to {8&y) come from regions where 
the distance between 7 and 7' is larger than £. 

We are now ready to calculate dephasing for those sys- 
tem paths shown in Fig. [21 As defined above, t\ and i 2 
are the two times at which the path 7' crosses itself. De- 
phasing acts on the loop formed by 7' and, as just argued, 
it acts once the distance between 7 and 7' is greater than 
£, i.e. in the time window from (ti+T^/2) to (t 2 — Tf/2), 
where T 5 (e) is given in Eq. (TT9]) , We average Eq. (|21b[) 
over the environment and use the central limit theorem 
to evaluate the action phase due to the coupling between 
system and environment, 

(exp[ic5$u]) onv = cxp [- ±(5$ u ) cnv 

= expHta-ti-Tf)/^], (27) 
where the dephasing rate is 

(28) 



h- 2 \-i y (u 2 



Given that (■ • • ) cnv is defined in Eq. (|22]) . we can sub- 
stitute Eq. directly into Eqs. ([2"T]). We thereby re- 
duce the problem to an integral over system paths which 
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is almost identical to the equivalent integral for U = 0. 
Assuming phase space ergodicity for the system, we get 
the probabilities 



(P sys (Ri, Y ;£i)) 



(P sys (Y,R 2 ;t-i 2 )) 



-ti/ru 



2na 



sys 

(t 2 -ti-T w /2)/T D 



2tt^ 



cos f e 



sys 

(t-*2-Tw/2)/T D 



(29a) 



(29b) 



(29c) 



2(Wl + W r )t d 

with il sys being the real space volume occupied by the 
system (the area of the cavity) . At this point the integral 
is the same as without dephasing, except that during the 
time (i 2 — £i — T^) , the inverse dwell time is replaced by 
(tq 1 +t^ 1 ). Thus when evaluating the (i 2 — ii)-mtegral, 
we get the extra prefactor exp [~T^{e) / t^] /(l + td/t^) 
compared with the equivalent integral result without de- 
phasing. Thus wc have 

e -T L (e)/T D -(T L (e)-T 5 (e))/T </l 

(F(Y ,e)) cx sine — . (30) 

Since (T L (e) - T s (e)) = r e with r c defined in Eq. (4}, 
the e dependence of the r^-term drops out. This means 
that (-F(Yo,e)) simply differs from its value without de- 
phasing by a constant factor, e -1 "^ 1 "* /(l + td/t^). Thus 
the integral over e in Eq. (|21a[) is identical to the one 
in the absence of the environment, and takes the form^ 
Re / °°dee 1+2 /( ATD ) exp[iE F e 2 /(\h)}, where we have as- 
sumed e <C 1. The substitution z = Epe 2 /(XH) imme- 
diately yields a dimensionless integral and an exponen- 
tial term, c~ Te / TD (neglecting as usual 0[l]-terms in the 
logarithm in t^). From this analysis, we find that the 
weak-localization correction is given by 

vl 



_,wl 



9o 



1 + Tt)/T0 



exp[-Tit/r^], 



(31) 



where <?q is the weak-localization correction at finite- Tp; 
in the absence of dephasing, 



5o wl 



■exp[-r^/T D ]- 



(32) 



We see that the dephasing of weak localization is not 
exponential with the Ehrcnfcst time, instead it is expo- 
nential with the AF-indcpendent scale given in Eq. ((4J. 
In all cases where £ is a classical scale (i.e. of similar mag- 
nitude to W, L rather than Af) we see that is much 
less than the Ehrenfest time, <C In such cases the 
exponential term in Eq. (|3ip is much less noticeable that 
the universal power-law suppression of weak-localization. 



E. Weak-localization for reflection and coherent 
backscattering 

We show explicitly that our semiclassical method is 
probability-conserving, and thus current-conserving, also 




Chaotic system 
environment not shown) 



Figure 3: (color online) A semiclassical contribution to co- 
herent backscattering for the system-environment model. It 
involves paths which return to close, but anti-parallel to them- 
selves at lead L. The two solid paths are paired (within W 
of each other) in the cross-hatched region. Here we show 
£ > eL, so the dephasing (dotted path segment) starts in the 
loop (Tj > 0). In the basis parallel and perpendicular to 
7 at injection the initial position and momentum of path 7 
at exit are r ± — (yo — y) cos 60, r || = (yo — y) sin #0 and 
Po± = Pf(0 - 6 ). 



in the presence of dephasing. We do this by calculat- 
ing the leading-order quantum corrections to reflection, 
showing that they enhance reflection by exactly the same 
amount that transmission is reduced. There are two 
leading-order off-diagonal corrections to reflection. The 
first one reduces the probability of reflection to arbitrary 
momenta (weak- localization for reflection), while the sec- 
ond one enhances the probability of reflection to the time- 
reversed of the injection path (coherent-backscattering) . 
The distinction between these two contributions is re- 
lated to the correlation between the path segments when 
they hit the leads. For coherent-backscattering contribu- 
tions, these segments are correlated (see Fig. [3]), but for 
weak-localization contributions, they are not. 

The derivation of the weak-localization for reflection 
r wl is straightforward and proceeds in the same way as 
the derivation for g wl given above, replacing the factor 
N K /(N R + N L ) by N L /(N R + N L ). We thus get, 



..wl 



1 + td/t^ 



exp[-r 5 /r 



(33) 



where r$ l = - ex-p[-T^/r D ] N[/(N L + N R ) 2 is the finite- 
Tg 1 correction in the absence of dephasing. 

We next calculate the contributions to coherent- 
backscattering, extending the treatment of Ref. [45| to 
account for the presence of dephasing. As before, the en- 
vironment is treated in the diagonal approximation. The 
coherent-backscattering contributions correspond to tra- 
jectories where legs escape together within Xw/2 of the 
encounter. Such a contribution is shown in Fig. [3J The 
correlation between the system paths at injection and 
exit induces an action difference <5<f>, 



sys 



SS c b s not given 
by the Richtcr-Sicbcr expression. It is convenient to write 
this action difference in terms of relative coordinates at 
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the lead (rather than at the encounter). The system ac- 
tion difference is then SS c f, s = — (po± + mAroj_)r j_ where 
the perpendicular difference in position and momentum 
are r 0J _ = {y - y)cos6 Q and p 01 _ = p F (0 - 6 ). As 
with weak-localization, we can identify three timescales, 
Tl,T^,T/, associated with the time for paths to spread 
to each of three length scales, L, W, £. However unlike 
for weak-localization we define these timescales as a time 
measured from the lead rather than from the encounter. 
Thus we have 

TKro±,Po±) - A- 1 rn[(mA£) 2 /|po± + mAr ±| 2 ] (34) 

with £ = {L, W,£}. Writing the integral over Yo as 
an integral over (roj_,po±) and using ppsin^odYo = 
dpoJ_dro±, the coherent-backscattering contribution is 



coupling, the environment part of Eq. (|12[) would then 
be 



J "CUV p p, 



1 + m- 



(38) 



cbs 



dpoxdrox ^r^i / F cbs (Yo) x _ (35) 
L pp sin uq 



After integrating out the environment in the same man- 
ner as for weak-localization, we get 

F chs (Y ) 

= (dY (' di(P sys (Y,Y o ,i))cxp[-(t-T^)/r ] 

N LPF sin O cxp [ - (Tl - T^/2)/rD - (T[ - T>)/r^ 



where \x is a number of order one, and A onv is the en- 
vironment level-spacing (for a two-dimensional environ- 
ment containing a single particle, A cnv ~ fi 2 /mL^ nv ). 
The first term above comes from the diagonal approx- 
imation used throughout this article, while the second 
term is a weak-localization correction. This correction 
becomes of order the system's weak-localization correc- 
tion on the timescale t ~ td , so there is a priori no reason 
to neglect it. 

The environment part of our calculation differs how- 
ever from the form-factor in that it corresponds to the 
time-evolution of the environment during the time it 
takes for a particle to be transported across the system. 
Therefore, the sum in Eq. (|38|) is not restricted to peri- 
odic orbits, and the unitarity of the environment's time- 
evolution imposes that fi = 0. Furthermore, unitarity 
(3glpst be preserved even in the presence of a finite-U, as 
long as there is no exchange of particles between system 
and environment. We thus conclude that we do not need 
to consider the weak-localization type corrections to the 
environment evolution because they cancel. 



n(N L + N R ) 



1 + td/t^ 



Now wc can proceed as for U = 0, pushing the mo- 
mentum integral's limits to infinity, and evaluating the 
tq±— integral over the range W , with the help of an Euler 
T-function. Wc finally obtain 



1 + td/t^ 



cxp[-r c /T ], 



(37) 



exp[-7£/T D ] N L /(N h +N R ), the finite- 
Tg 1 coherent backscattering contribution in the absence 



in terms of r{j s 



of dephasing. Hence r + r = —g for all values 
of and Ttf>, and our approach is probability- and thus 
current-conserving. 



F. Weak-localization corrections in the 
environment 

So far, we have only considered cases where we make 
a diagonal approximation for the environment. On the 
face of it this seems a little unreasonable. For instance, if 
the system and environment are of similar sizes then one 
would expect that diagonal contributions for the system 
and weak-localization for the environment would be as 
important as the contributions calculated above. 

Since the environment is a closed cavity, one would 
naively think that the weak-localization contribution for 
the environment should be calculated in a similar manner 
to the form-factor in Refs. I7lll74ll. In the absence of 



Regime of validity of the semiclassical 
calculation 



Throughout this article we assumed that the system- 
environment coupling is weak enough not to modify the 
classical paths in the system. Formally, this assump- 
tion can be rigorously justified by invoking theorems on 
structural stability^. However, care should be taken in 
extrapolating our results to the limit £ — ► 0, since the 
force on the particle is the gradient of the interaction po- 
tential, ~ U/£. We therefore estimate the minimum £ for 
which we can legitimately assume that classical system 
paths are left unchanged by the system-environment cou- 
pling. This will give the bound on the regime of validity 
of our approach. 

To see significant dephasing we need ~ td, so we 
cannot take the interaction strength to zero, instead we 
require that (U 2 ) ~ X cnv h 2 /TD, see Eq. ([28]). This in- 
duces a typical force ~ (tl 2 ) 1 / 2 /^ ~ (hj '£)(A cnv /td) 1 / 2 on 
the particle. To see if this noisy force significantly modi- 
fies the paths in the vicinity of the encounter, we compare 
it with the relative force of the chaotic system Hamilto- 
nian on the particle at the encounter. Since the perpen- 
dicular extension of the encounter is Sr± ~ (LAf) 1 ^ 2 , 
and the duration of the encounter is of order the Lya- 
punov time ~ A -1 , the system force goes like mX 2 5r± ~ 
toA 2 (LAf) 1 ^ 2 - Estimating A -1 ~ v-p/L as is typical of 
chaotic billiards, the ratio of the noisy force to the sys- 
tem force becomes [X ciiv LXf / '(£ 2 X 2 td)]~ 1 ^ 2 . Thus one 
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Chaotic system 
(environment not shown) 

Figure 4: (color online) Dephasing of weak-localization when 
^<Le~ (LAf) 1,/2 , and hence T e (e) is negative, see Eq. (fl9)l . 
The dephasing starts and ends in the "legs" rather than the 
loop. In the language of disordered systems this means the 
dephasing affects the diffusons as well as the cooperon. 



can ignore the modifications of the classical paths due to 
the coupling to the environment, as long as 



t » (A F £) 



1/2 



Aonv/A 



1/2 



(39) 



We see that £ can easily be less than the typical encounter 
size Sr± ~ (LAp) 1 / 2 (remember that Atd ^ 1 is always 
assumed). Thus our method is not only applicable for 
£ up to the system size, where dephasing happens only 
in the loop. It is also applicable for £ smaller than the 
encounter size, in which case the time T^(e) is negative, 
and dephasing occurs in part of the legs as well as the 
whole of the loop, see Fig. |U 

Finally we caution the reader that the whole semi- 
classical method used in this article relies on the lead 
width being greater than the encounter size, this requires 
that Atd <C (L/Af) 1 / 2 , thus we cannot access the regime 
£ ~ Af, which is dominated by stochastic diffraction at 
the leads. 



H. Shot noise in the presence of an environment 

When the temperature of the electrons in the leads 
coupled to a chaotic system is taken to zero, there is no 
thermal noise in the current through the device. However 
there is still the intrinsically quantum noise which orig- 
inates from the wave-like nature of the electrons. This 
zero-temperature noise is known as shot noise^. In the 
absence of dephasing, shot noise has been well-studied 
using RMT— , quasi-classical field theory^, and semi- 
classical method o 49 i 50 ' 76 . 



It is generally argued that the shot noise is unaffected 
by the presence of an environment which causes dephas- 
ing but not heating of the electrons - the regime of 
phase-breaking of Ref. (5?| . This belief is founded on the 
fact that (i) the dephasing-lead model gives a dephasing- 
independent shot noise, (ii) kinetic equations - in which 
interference effects are ignored - give the same shot noise 
as full quantum calculations. Here we show explicitly 
that, under the assumption that the system-environment 
coupling does not heat up the current-carrying electrons, 
indeed, the coupling to the environment does not affect 
shot noise. 

We start with the formula for the zero-frequency shot 
noise power through a system coupled to an environment. 
This formula is derived in Appendix IA 21 and is given as 
Srr(0) in Eq. (|A17|) . We use Eq. (JTTJ) to write each ma- 
trix element as sums over classical paths. This gives us a 
sum over eight paths - 4 system paths and 4 environment 
paths - as sketched in Fig. [5l The system paths are as 
follows 



71 from 2/01 on lead L to 2/1 on lead R, 

72 from ?/03 on lead R to 2/1 on lead R, 

73 from ?/03 on lead R to 2/3 on lead R, 

74 from 2/01 on lead L to 2/3 on lead R. 



The sums over lead modes and the trace over the environ- 
ment density matrix are performed in the same manner 
as for the conductance, [see above Eq. (|12p]. which results 



e*V 
WE? j 



/ dy 01 / d2/03d2/id2/3 V] AjysC 1 * 8 



/ dqoidq 03 dqidq3 ^ A 



71.---74 

(*c„v+* lt ) (4Q\ 



n,---r4 



Here, A s 



A 1 \A 1 2A 1 j,A 1 ^ 



4>s 



(S. 



71 



S 



-2 



5 7 3 — S 1 A c )/h and we absorbed all Maslov indices into 
the actions. Similarly, A env = A^\ Ar2^4 F 3^-r4i ^nv = 
(Sri - S T 2 + Sr3 - Sn)/h, and $ u = (§ 7 i,n - § 7 2,r2 + 
§73, T3 — §74,r4)/^- As argued above, in the regime of 
pure dephasing, <& S ys, ^onv and $u are uncorrelated. We 
thus first pair up the system paths to minimize $ S ysi this 
pairing is the same as it would be in the absence of the 
environment (compare Fig. [5] to Fig. 1 of Ref. [E(|). We 
see from the construction of Eq. (|A17[) , that all paths 
reach the encounter at the same time££, t\ . 

Now we make the crucial observation that, for any 
given set of system paths, we have 7I ~ 72 for times 
greater than t[. Thus for times greater than t[ we can 
write the sum over Tl, T2 in the second line of Eq. (f¥D|) 
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Figure 5: (color online) Sketch of typical trajectories which 
contribute to shot noise in the presence of an environment. In 
both (a) and (b) the system (environment) paths are sketched 
above (below) the time axis. In (a) we show a contribution 
which will survive system averaging, because the system paths 
are paired almost everywhere with an encounter at time t[. 
There is no constraint on the environment paths, as yet. For 
simplicity we show only system paths 7I and 73, with path 
72 and 74 being the same as 7I and 73 except that they cross 
at the encounter. Depending on the choice of t'\ with respect 
to the other timescales, ti,t2,i3, this pair of system paths 
could represent any of the contributions in Fig. 1 of Ref. [Hoi ]. 
In (b) additional constraints are imposed on the Pi's, after 
one has integrated over all initial and final positions of the 
environment. This integration removes all contributions from 
environment propagation when the system paths are paired. 



I dqi J2 AriA r2 c i{Sri - Sr2+$ -> 1 ' ri - $ ~< 2 - r2)/h 
J n,r2 

- / dqi J2 AriAr 2 e [{Sri - Sr2+$ -> 1 - ri ~ s -' 1 - r2)/h 
J n,r2 

= J dqj [X' env (q 1 ,t 1 ;^(,t[)]*X' env (q l ,t 1] ^ 1 ,t' 1 ) (41) 



where we set 72 = 7I to get the second line. We de- 
fine 3Cg nv (qi, t%; q^,ti) as the propagator for the environ- 
ment evolving under an effective time-dependent poten- 
tial, V'(q, t), which is the sum of U(r 7 i(i), q) and the po- 
tential term in H cnv . Since both propagators in Eq. (|41[) 
evolve under the same potential (because 7I = 72), we 
can use the basic property of propagators^ that 

J dqi [3C^ v (qi,t 1 ;q / ,ti)]*3C^ v (q 1 ,ti;q ,/ ,ti) 

= S(q" - q') (42) 

to integrate out these propagators for times greater than 
t\. In their place we have a constraint that path Tl and 



T2 must be the same at time t[. This is sketched in 
Fig. the paths IT and T2 to the right of the encounter 
in Fig. [5^ are replaced in Fig. [5b by the constraint that 
paths n and T2 meet at time t[. Note that we cannot use 
Eq. (|42|) to integrate out paths Tl and T2 for arbitrary 
times before t[ because the system paths 7I and 72 are 
then different enough that the potential V'(q,t) will be 
different for the two propagators in Eq. (|41[) . However 
we can use the same argument to integrate out the pair 
r3-r4 after time t\, and to integrate out the pairs ri-r4 
and r2-r3 before the time t\. After all these pairs are 
replaced by 8- functions, we get the situation shown in 

Fig.rjb. 

Focusing on £ much greater than the encounter size, 
the above method can be used to integrate out the envi- 
ronment paths for all t > t\ and all t < t[. This leaves a 
single point (the environment state at t = t^) to integrate 
over. Doing this we see that Eq. (|40|) reduces to 



(43) 



S = tt-tt^ / dyoi / d2/ 03 dyidy 3 ^ -W,«-" : ' > 



This is identical to the shot noise formula in the absence 
of an environment. Thus we have completely removed 
the environment from the problem without affecting the 
shot noise of the system at all. To calculate the shot 
noise now, one simply needs to follow the derivation for 
a system without an environment in Ref. (50j . The result 
is most conveniently written in terms of the Fano-factor, 
F, which is the ratio of the shot noise to the Poissonian 
noise 2e(I), where (I) = 2e 2 g D V/h is the average current. 
One gets 



S/2e(I) 



N h N R 



(N L + N R )- 



■exp[-r° p /r D ]. (44) 



This is of course independent of the coupling to the en- 
vironment. 

As a final comment, we note that above we kept 
only the leading 0[N] term in the shot noise. There is 
a hierarchy of weak-localization-like corrections 0[7V a ], 
a = 0, — 1, . . . to this resuh— , which are suppressed by 
dephasing in much the same way as the weak-localization 
correction to conductance. Thus for Tg P < td, we can 
expect the environment to cause a cross-over from the 
result in Ref. [zi] to the result in Eq. gj) with r° p = 0. 
Hence in the classical limit of wide leads (A^l,r ^> 1) 
the environment's effect is negligible, however for nar- 
row leads (-/Vl,r ~ 1) the environment's effect may be 
significant. 



III. CLASSICAL NOISE 

We have shown that, to capture the effect of dephas- 
ing on weak localization, it is sufficient to treat the en- 
vironment at the level of the diagonal approximation. 
We thus observe that, in the semiclassical limit of short 
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wavelength, Xp/L env — > 0, a quantum chaotic environ- 
ment has the same dephasing effect on weak-localization 
as the equivalent classical chaotic environment. Because 
correlations typically decay exponentially fast in classical 
hyperbolic systems, this makes the effect of this classical 
environment very similar to a classical noise field with 
a suitably chosen spatial and temporal correlation func- 
tion. In this Section we show that the conclusions that 
we draw for a quantum environment can also be drawn 
for a classical noise field. One common experimental 
example of such a field is microwave radiation, applied 
to the chaotic dot either by accident or on purpose^!. 
A second example, is the common theoretical treatment 
of electron-electron interactions as a source of classical 
(Johnson-Nyquist) noised. 

We add a new term to the system Hamiltonian of the 
form, V no i se (t). We assume this term is weak enough that 
it does not affect the classical paths, but strong enough to 
modify the phase acquired along such paths. The phase 
difference for a pair of paths contributing to the conduc- 



tance is ($ 
and 



i i 



sys 



dt[V(r Y ;t)-V{r y ;t)]/h. 



(45) 



We now assume that the noise is Gaussian distributed 
with 

(Vr(ry(t);t)V(r 7 (tV)) 

= (V 2 ) K noise (\r 2 - n|/0 J„oi S c(A no i S c|i2 - til). (46) 

Here, K no i se (x) gives the form of the spatial decay of the 
correlation function (on a scale £), and J no ise(t) gives the 
form of the temporal decay of the correlation function (on 
a scale which we call A~ oise to make the analogy with the 
notation in Section III C[) . We can now follow the deriva- 
tion in Section III C[) by replacing U(r 7 (i), qr(i)) with 
V(r 1 \t) throughout. We assume that the correlations in 
time are short enough to be treated as white-noise-like, 
<^noise(aO oc 5(x), and that the spatial correlations de- 
cay fast enough that we can justify the approximation 
in Eq. (|26p . This directly leads to the same result for 
weak- localization as in Eq. (|3ip . where now 



2\-l 

noise 



{v 2 



(47) 



A. Noise due to electron-electron interactions in a 
2-dimensional ballistic system. 

Here we consider the noise generated by electron- 
electron interactions in a ballistic chaotic system. In such 
a system, dephasing is caused by noise with momenta (8p- 
vectors) larger than the inverse system size, L . Thus 
the dephasing processes are the same as those for ballis- 
tic motion in disordered systems (Sp- vectors greater than 
inverse mean free path). Such processes were first stud- 
ied in Ref. [l3j for spinless electrons, while more recently 



Ref. [l5j explored the full crossover from ballistic to dif- 
fusive motion for electrons with spin. The effect of a 
finite Ehrenfest time on such dephasing was considered 
in Ref. [24j, which found that the dephasing rate in the 
vicinity of the encounter has a logarithmic dependence on 
the perpendicular distance between paths. This led them 
to observe that the electron-electron interaction in a 2- 
dimensional ballistic system dephases weak-localization 
exponentially with the Ehrenfest time. We repeat their 
derivation here and show that 



wl 

9a 



f + td/t^ 



exp [- (t£ + ±t Lt ) / T<t> ] , (48) 



), where <& sys is given in Eq. (flSa) Rcfs. 



where tx t is given by Eq. (|3ip withf equaling a thermal 
length scale Lt = hv-p/k^T. Ref. [2J] neglected the tl t - 
term in the exponent since it is often small. 

In our qualitative derivation of this result, we treat 
the electron-electron interaction as classical-noise, rather 
than using the perturbative field theory approach in 
T3lfl5l .24]. Our approach is similar in spirit to those 
for diffusive systems^. The screened electron-electron 
interaction gives a correlation function of the formic 



(V? 



L P ■ Sp) 



\Sp\ 



x J{uj), 



(49) 



where we assume |<5p| -C |p|, so the energy difference 
between a particle with momentum (p + Sp) and p is 
p • 5p/m. The factor of l/|<5p| comes from the imagi- 
nary part of the screened Coulomb interaction, is due to 
the polarization bubble and corresponds to the fluctua- 
tions of the electron sea at momentum and energy (Sp, u>). 
The 5-function ensures that energy and momentum are 
conserved in the interaction between the system and the 
environment. The function J(to) gives the weight of en- 
vironment modes excited at energy u>. At temperature T 
it is typically of the form J(lo) ~ [ sinh(o>/fcBT)] ■ It 
is convenient to write Sp in terms of components parallel 
and perpendicular to the relevant classical path, i.e. par- 
allel/perpendicular to p. Then, 



vpSpu 



[5p\ + Sp\Y^ 



x J(oj), 



(50) 



and we have 

(V(r Y (t');t')V(r 7 (t);t)) 

= J d d Sp J d^psrd'^+uit'^/n (v£ j(51a) 

oc / duj J(ui) exp [i2u(t' -t)/h] 

cxjp[iSpj_5r±/h] 



xRe 



l/L dSp ^{spi + (u/v F n/i 



(51b) 



Here, Sr(t',t) = (ry(t') — r 7 (i)). To get the second line, 
we wrote Sr(t',t) = (5r\\,5r±) in the basis parallel and 
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perpendicular to p, we then inserted Eq. (|50|) . evalu- 
ated the (5p||-integral and noted that Sr^ = vp(t' — t). 
Eq. ([5"Tj) expresses the correlator of the Coulomb inter- 
action along classical trajectories in the form convenient 
for our semiclassical approach, Eq. (|46|) . At first sight 
the form of the interaction in Eqs. ()49|50|) does not ap- 
pear to correspond to a classical noise-field, since it is 
a function of the momentum, p 7 (t), of system paths. 
The correlator must however be evaluated on weak lo- 
calization loops, in which case one can use the fact that 
r||y(i) = ?"|| 7 (i) = v-pt throughout the encounter region 
to perform the Fourier transform. Thus the interaction 
term in Eq. (|51b[) is equivalent to a classical-noise field 
which is a single function of <5r and t for all p 7 . This 
function must simply be chosen such that the integral 
over Sp» reduces to Eq. (|51b[) . 

The time-dependent part of the correlator is given by 
the w-integral in Eq. (|51b|) . We assume that the tem- 
perature is high enough, k^T > hvp/L, that the corre- 
lation time becomes shorter than the time of flight L/vp 
through the cavity. Accordingly, we treat the noise as 
^-correlated in time, and set t' = t from now on. 

We next investigate the properties of the real part of 
the <5p^-integral in Eq. (|51b[) . giving the spatial depen- 
dence of the correlator. We write it as 



G(5r±_) = Re 



p F Sr ± /h 



d.c c 1: 



5r±/L [* 2 + (SrjL u y] V2 



(52) 



where = hvp /lu is the distance a system particle will 
travel on the timescale that the w-cnergy component of 
the noise fluctuates. For the ballistic model of the e- 
e interactions to be valid we need that L^ <C L, but 
we assume that ^> Xp. We can easily evaluate this 
integral in the following regimes^, 



ln[Lx/AF] — ln[w//cBT]. We can neglect the second term 
as it is much smaller than the first, and thereby replace 
L u with Lt in the above formulae. In this way we re- 
produce the result in Ref. 24[, that the dcphasing rate 



in the vicinity of the encounter is 



while f, 1 (Sr± > Lt) 



__ 1 lnjJrj_/A F j 
* ln[L T /A F ] 



(55) 



This is rather differ- 



ent from the systems considered elsewhere in this arti- 
cle where the dcphasing rate is approximately zero for 
5r± < £ and approximately constant for 5r± > £. 

We now calculate the effect of such a (Jr^-dcpcndent 
dephasing rate. We note that close to the encounter, 
5r± = ^eLc Xr where r is the time measured from the 
encounter. We split the dephasing into two contribu- 
tions. The first contribution is where Sr± > Lt (here 
dephasing is time-independent, at the rate t^ 1 ), the sec- 
ond is where the paths have Ap < Sr± < Lt (here 
dephasing is time-dependent). The boundary between 
the two contributions is at r = Tr(e)/2, where we de- 
fine T T (e) = A" 1 ln[L^/(ie) 2 ]. The lower bound on 
the second contribution (5r± = Af) is at time r = 
Xr(e)/2 — X^ 1 \u[Lt/ Xp]. Then the exponent induced 
by the dephasing is 



U 



Me) 



TV, 

t 2 -h 



T T /2 



dr 



Me) 



TV, 



Tt/2-A- 1 1ii[L t /A f ] t 4> 

X- 1 ln[L T /A F ] 

T 4> 



(56) 



ln[L w /Ap], for <5rj_ <c Af, 

G(Sr ± ) ~ { ln[L u /Sr x ], for A F < Sr x < L u , (53) 
0, for L u <C 8r± ~ L, 

where we neglected all 0[l]-terms (for example the re- 
sult for Sr ± > L u is actually 0[L u /6rx] < 0[1]). The 
crossover between these regimes is smooth. We thus con- 
clude that K no i se (x), as defined by Eq. (|46|) . becomes 



A, 



,(6r/L„) = 



ln[£^/frj 
ln[L u /X F 



(54) 



This function does not decay fast as Sr± grows, and thus 
it cannot be treated in the manner we do elsewhere in 
this article, i.e. we cannot write an equation analogous 
to Eq. (|26|). Instead we note that the dephasing rate in 
the vicinity of the encounter (where Af <C 5rj_ <C L u ) 
goes like [G(0) - G(A F < Sr ± < L u )] = ln[<5r ± /A F ]. 
While the dephasing rate in the loop, t7 , goes like 
[G(0) - G(Sr ± ~ LJ)] ~ ln[L w /A F ]. Now we note that 
the integral over u> is dominated by us ~ k^,T <C Ep, thus 
we can define Lt = hvp/kp$T and write ln[L w /Ap] = 



where to evaluate the integral we defined r' = r — 
Xr(e)/2. The first term (which comes from the dephas- 
ing in the loop) alone would give no exponential term 
in the dephasing. The integral of that term over ti — t\ 
gives it a form Tl(c) — Tt(c) = tx t giving an exponent 
like in Eq. (f3"Tj) with £ = Lt- However we also have 
the second term which gives dephasing in the vicinity of 
the encounter, we can write it in terms of an Ehrenfest 
time using A -1 h^Lx/Ap] = — |t£, t . Summing the 
two terms we find that the exponential term in the de- 
phasing goes like f/r^, with f = rpj + |tx t . This is 
the result which we gave in Eq. (|48[) and was found in 
Ref. [24| (neglecting the 7x T -term). 

Because Af is the scale of Friedel oscillations, one 
might have expected that electron-electron interactions 
give a noise with a correlation length £ ~ Af, which 
would lead to a suppression cx exp[— 2t£c /t^\ of weak lo- 
calization, instead of exp[— r| /r^]. The factor of 2 differ- 
ence between the correct result and this naive argument, 
is due to the fact that all scales (i.e. all values of Sp) 
contribute to the noise induced by the electron-electron 
interactions. 
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3 dimensional metal gates 




Figure 6: Sketch of a typical chaotic quantum dot in a two- 
dimensional electron gas (2deg). The dot is defined by metal- 
lic gates which are biased to deplete the 2deg everywhere ex- 
cept in the regions defining the leads and the dot. These 
gates are a distance D above the 2deg (with D S> A| atc ). At 
finite temperature, the electrons at the surface of these metal- 
lic gates will fluctuate, leading to noise which will be felt by 
the electrons in the chaotic quantum dot. 



trap at the edge of a gate - in which case their noise 
field will be similar to that in Ref. This modifies 

the integrand of the g-integral, however the upper cut- 
off will still be h/D, therefore dephasing will again be 
Tg-independent. 

In general, dephasing is due to a combination of e-e in- 
teractions within the system and the Coulomb coupling 
between the system and external charge distributions in 
gates and other reservoirs of charges. Which of these 
sources of dephasing dominates is determined by micro- 
scopic details which we do not discuss here, in particu- 
lar the temperatures and mean-free-paths of both system 
and gates>^. 



IV. DEPHASING LEAD MODEL 



B. Noise due to the coupling to the electrostatic 
environment 

Our aim here is to show that electron-electron interac- 
tions do not automatically lead to dephasing which is an 
exponential function of the Ehrcnfcst time. It only hap- 
pens when the g-integral is divergent at its upper limit, 
with an upper cut-off of order pp- If the (/-integral is 
cut-off by some other length scale, then the dephasing 
of weak-localization will be independent of the Ehrcnfcst 
time. 

The example we consider is noise in a two-dimensional 
system due to electron-electron interactions between the 
system and the gates. A typical experimental set-up is 
sketched in Fig. [6] The gates are bulk metal, the elec- 
trons in the system will feel fluctuations of electrons at 
the surface of the gates. One can expect that these fluc- 
tuations at the surface of the gate are sufficiently well- 
confined to two-dimensions (by screening in the bulk of 
the gate), that they will cause a noise field in the system 
of a form similar to Eq. (|49p . However the fact that the 
distance between the gates and the chaotic system is D 
means that the natural upper cut-off on the q-integral will 
be h/D and not p-p. Assuming fl«L, wc replace pp 
by h/D throughout the derivation in Section fill A[ and 
find that 



wl 
i/gatc c 
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1 + TD / T0 



cxp[-t c /t ], 



(57) 



with £ = (LtD) 1 / 2 and hence t c = A" 1 ]n[L 2 /(L T D)]. 
The dephasing rate here is t^ 1 oc D~ x . In systems in 
which dephasing is dominated by thermal system-gate 
interactions, wc therefore expect a dephasing that is te- 
independent and algebraic in for D ~ L. 

In real experiments the gates are typically much more 
disordered than the chaotic system, thus we can easily 
have a situation in which the thermally excited modes 
in the gates are diffusive - or even localized by a charge 



In its simplest formulation the dephasing lead model 
consists of adding a fictitious lead 3 to the cavity. This 
is illustrated in Fig. [7J Contrary to the two real leads 
L,R, the potential voltage on lead 3 is tuned such that 
the net current through it is zero. Every electron that 
leaves through lead 3 is replaced by one with an unrelated 
phase, leading to a loss of phase information without loss 
of current. 

In this situation the conductance from L to R is given 
byI6 



9 



RL 



?3l + Tm ' 



(58) 



where T nrn is the conductance from lead m to lead n 
in the absence of a voltage on lead 3. We separate the 
Drude and weak- localization parts of T nm , 



rji npD 

-Lnm — 1 nm 



6T n 



(59) 



where the Drude contribution, T„ m , is ©[A^t] and the 
weak- localization contribution, ST nm , is 0[A^] and = 
Al + Ar + Na is the total number of channels in this 
three terminal geometry. We expand g for large A^r 
and collect all ©[Ay-terms (Drude contributions) and 
all 0[AT.j] -terms (weak-localization contributions) to get 
q = q° + o wl with 



rpD 

RL 



J 3L "r 1 



3R 



RL 



(T^ST 3h + (T^5T 3R 



(T. 



3R 



3L7 



(60a) 



(60b) 



These equations form the basis of our scmiclassical 
derivation of weak localization in the dephasing lead 
model. We first consider the case of a dephasing lead 
perfectly coupled to the cavity, and then move on to con- 
sider a dephasing lead with a tunnel barrier of trans- 
parency p < 1. We finally discuss multiple dephasing 
leads. 
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adR 



lead 3 (dephasing) 

Figure 7: Schematic of the dephasing-lead model. The system 
is a open quantum dot with an extra lead (lead 3), whose 
voltage is chosen to make zero current flow (on average) in 
that lead. This extra lead thereby causes dephasing without 
loss of particles. 



A. Dephasing lead without tunnel-barrier 

1. Weak-localization 

The Drude conductance and weak-localization correc- 
tion from lead m to lead n in a three-lead cavity are 



T 



N n N m 

N n Nr, 



■cxp[ 



/td] 



(61) 
(62) 



where f^ 1 = (t q L)^ 1 (W^ + W R + W 3 ), in terms of r , 
the time of flight across the system. We substitute these 
results into Eq. (|60a|) and Eq. (|60b|) , and write the answer 
in terms of the dwell time in the two lead (L and R) 
geometry, td, and the dephasing rate 
define as the decay rate to lead 3, 



Tj, 1 . which we 



= {t q L)- 1 {W^ + Wk), 



(63) 
(64) 



We have f D 1 



t, , and from this we find the 



Drude conductance and weak-localization correction, 



, wl 



,.,wl 



1 + Tu/Ttj, 



expf-TE 1 /^]. 



(65) 
(66) 



Here and g^ 1 are the results for a two lead cavity in 
the absence of dephasing, Eqs. p7|) and (|32|) 

The weak-localization correction with a dephasing lead 
has a similar structure to that with a real environment. 
However here the time scale involved in the additional 
exponential suppression contains no independent param- 
eter analogous to £. We could have expected that the 
width of the dephasing lead would play a role similar to 
£. However this turns out not to be the case, instead the 



Fermi wavelength appears in place of £, so the time scale 
in the additional exponential suppression is the Ehrenfest 



time, r, 



E 



2. Universal conductance fluctuations 

In the absence of tunnel barrier, we can go further 
and calculate conductance fluctuations at almost no ex- 
tra cost. From Eqs. (|58[) and ([59]) . we get the following 
expression for the variance of the conductance to order 
0[A°] 



var g 



var T RL 

(T3D) 4 var T 3L + (T 3 D L ) 4 var T 3R 



^3 D l) 4 



(^3R ^3LJ 

(T 3 D R T 3 D J 2 covar (T 3L ,T 3R ) 



(r 3 u L ) 2 covar(r RL ,T 3R ,) 



(^3R "T J 3L^ 

(T 3 D R ) 2 covar(T RL ,r 3L ) 



3R 



(67) 



In the universal regime, Ref. [T(| gives, to order 0[N%] 
in the presence of time-reversal symmetry, 



varTL 



covar (T y ,T ifc ) 



NfN 3 N k 



(68) 
(69) 



Ref.[5i 

te/td- 



showed that Eq. |68|) remains valid even at finite 
Inspection of their calculation for varT^ shows 
that the same conclusion also applies to covar (Ty, Tjfc), 
and thus Eq. (|6^|) still holds independently of te/td. To- 
gether with Eq. (f6"3")) . this straightforwardly leads to 



vai g 



N*N* 



1 



(iV R + iV L ) 4 (l + r D /r ) 5 



(70) 



In the universal regime, this result was previously derived 
in Ref. [ll|. We thus conclude that, in the dephasing lead 
model without tunnel-barrier, conductance fluctuations 
exhibit the universal behavior of Eq. (|70|) . Below, we 
confirm this result numerically. 



B. Dephasing lead with tunnel-barrier 

Putting a tunnel-barrier on the dephasing lead 3 is at- 
tractive because one can avoid the local character of the 
dephasing lead model by considering a wide third lead 
with an almost opaque barrier—. Additionally, this is 
the model studied numerically in Ref. [23[ in the con- 
text of conductance fluctuations. Weak-localization and 
shot noise in this model have been considered within the 
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trajectory-based semiclassical approach in Ref. [56j, and set of equations 
here we only mention the main results. 

According to Ref. [f56| , when all leads are connected to Ir 
the cavity via tunnel-barriers, with the barrier on lead m = 1 

having transparency p m € [0, 1], the Drude conductance 



-T sub V 



+ T L y L 



(74a) 
(74b) 



between lead m and n, 
correction, ST nm7 are 



T„ m , and the weak- localization 



nm 



N n N„ 



(71) 



PnPr, 



Pn + Pm - j7 



x exp [-r° p /r D2 - (r$ - t° p )/t Di ] . (72) 

Here, r^ 1 = (j^L)- 1 J2 n p n W n and = 
(tqL)^ 1 J2 n PnQ - Pn)W n arc the single path and 
the paired paths survival times respectively, W n is the 
width of lead n, N = J2k PkNk and N = J2k Pk^k- 

Now we assume that pw = PR = 1 so only the dephas- 
ing lead has a tunnel barrier. Substituting the Drude 
and weak-localization contributions into Eq. (|60b[) we 
find that 



. wl 
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5o wl 



1 + Tux/To 



exp [-(1 



p) t e/ t <I> 



fA*] ■ (73) 



The argument in the exponent in Eq. (|73[) has a simple 
physical meaning. It is the probability that a path sur- 
vives throughout the paired-region (t^P /2 on cither side 
of the encounter) without escaping into lead 3, multiplied 
by the probability to survive the extra time (r^ 1 — t^ p ) 
unpaired without escaping into lead 3 (to close a loop of 
length Tg 1 ). The first probability is exp[— (2 — p)t£F /t^] 
while the second is exp[— (t^ 1 — tP-P)/^]. 

We note that if we consider a nearly opaque barrier, 
the relevant time scale involved in the exponent is t^ 1 + 
Tg P ~ 2tjq . Thus by tuning the opacity of the barrier, we 
can vary the exponential contribution to dephasing from 
cxp[— Tg 1 /t^] to exp[— 2Tg /t^], but we cannot remove the 
exponent. In particular, we cannot mimic dephasing due 
to a real environment with £ ~ L, Eq. (f3Tj) . since it has 
only a power-law dephasing. 

There is to date no theory for conductance fluctua- 
tions at finite Ehrenfest time in the presence of tunnel 
barriers, and constructing such a theory would require a 
formidable theoretical endeavor. Numerically, Ref. [23| 
observed that, in contrast to the dephasing lead model 
without tunnel-barrier, conductance fluctuations in the 
presence of a dephasing lead with a tunnel barrier are 
exponentially damped oc exp[— te/t^]. 



where the superscript-T indicates the transpose. The 
column-vectors I, V have an ith element given by the 
current or voltage (respectively) for the dephasing lead 
i € {3,n + 2}. The column- vectors Tl and Tr have an 
ith element given by Tn and Tm, respectively. Finally 
T su b has matrix elements given by 



(75) 
(76) 



where again i,j G {3, n + 2}. Substituting V from 
Eq. (|74bj) into Eq. (|74aj) and using /r = <j>Vl gives us 
the conductance from L to R as 
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r, 



LR 



L i S ub i R 



(77) 



Thus finding g requires the inversion of T su b- This is 
cumbersome, so instead we present a simple argument to 
extract only the information we are interested in - the 
nature of the exponential in the dephasing. 

We argue that whatever the formula for conductance 
for n dephasing leads is, we can expand it in powers of N 
and collect the 0[7V°]-terms to get a formula for weak- 
localization of the form 



71 + 2 



ST LK + J2 Aj5T Lj + B 3 5T ]K + 



3=3 



n+2 

E 



CijSTij.{79) 



Here, the sum is over all dephasing leads. To get the pref- 
actors Aj , Bj , Cy we would have solved the full problem 
by inverting T su b, however we can already see that, to 
leading order, they are combinations of Drude conduc- 
tances and thus independent of the Ehrenfest time. In 
contrast all the weak-localization contributions contain 
an exponential of the form [td 1 and td 2 are defined be- 
low Eq. d?!])] 



exp[ 



op 



Ad 2 + (te -t e p )/t Di 



(79) 



Thus defining r^" as the rate of escaping into any of the 
dephasing leads, so rT" 1 = (t L) -1 Sj=3 PjWj, we see 
that <7 wl decays with an exponential 



g w oc exp[-(l - /5)te P /t 



(80) 



C. Multiple dephasing leads 

The n probe dephasing model consists of adding n fic- 
titious leads to the cavity (labeled {3, • • • ,n+ 2}) in ad- 
dition to lead L, R. The voltage on each supplementary 
lead is tuned so that the current it carries is zero. With- 
out loss of generality we defined Vr = 0, then we get the 



where we define p such that pr^ 1 = (r i) _1 J2j PjWj ■ 
We have just shown that multiple dephasing leads 
cause an exponential suppression of the weak-localization 
which is qualitatively similar to that caused by a sin- 
gle dephasing lead. The exponent is proportional to the 
Ehrenfest time, and contains no independent parameter 
analogous to £. 
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NUMERICAL SIMULATIONS 



A. Open kicked rotators 

Because of the slow, logarithmic increase of te with 
the size M of the Hubert space, the ergodic semiclassical 
regime te > td, Atd 3> 1 is unattainable by standard 



numerical methods. We therefore follow Refs. [23II44II4S 
l5lll52ll8"ol | and consider the open kicked rotator model. 



1. Dephasing lead model 

The system is described by the time-dependent Hamil- 
tonian 



H = 



(p - Po? 



K cos{x — xp) 8{t — nrp). (81) 



The kicking strength K drives the dynamics from inte- 
grable (K = 0) to fully chaotic (K > 7). The Lyapunov 
exponent in the classical version of the kicked rotator is 
given by At « \n(K/2). In most quantum simulations, 
however, one observes an effective Lyapunov exponent 
A e ff instead that is systematically smaller than A by as 
much as 30%21. Two parameters, po and Xq, are intro- 
duced to break two parities and drive the crossover from 
the /3 = 1 to the (3 = 2 universality class^I, corresponding 
to breaking the time reversal symmetry. 

We quantize the Hamiltonian of Eq. (f8~Tj) on the torus 
x,p G [0, 2w], by discretizing the momentum coordinates 
as pi = 2ttI/M, 1 = 1,... M. A quantum representation 
of the Hamiltonian of Eq. (I81|) is provided by the unitary 
M x M Floquet operator U, which gives the time evolu- 
tion for one iteration of the map. For our specific choice 
of the kicked rotator, the Floquet operator has matrix 
elements 



Vi. 



M -l/2 e -(m/2M)[(l-l„) 2 + (l'-lo) 2 ] 



v 2irim(l-l')/M-(iMK/2ir) cos[27r(m-m )/-M1 

x e e 1 



with Z = p M/2n and mo = xqM/2tt. The Hilbert space 
size M is given by the ratio L/X-p of the linear system 
size to the Fermi wavelength. 

To investigate transport, we open the system by defin- 
ing n + 2 > 2 contacts to leads via absorbing phase-space 
strips [k - N t /2 1 k + Ni/2 - 1], i = 1,2 . . .n + 2. With 
iVx = J2i Ni, we construct a Nt x Nt scattering matrix 
from the Floquet operator U 



s(e) = V3~¥y<-y e- iE o-uV3^~ ¥y uy f . (83) 

The Nt x M projection matrix D 5 , which describes the 
coupling to the leads, has matrix elements 



O 5 



JlSi m iHGUi^}, 
otherwise, 



(84) 



where \-fi\ 2 = pi G [0,1] gives the transparency of the 
contact to the external channel I, and |L {I 1 } denotes 
the ensemble of cavity modes coupled to the external 
leads. Below we focus on perfectly transparent contacts, 
Pi = li = 1, VL 

The conductance (|60a[) in the dephasing lead model is 
obtained from Eq. (|83")) with n + 2 = 3 leads. The trans- 
port leads carry N = N^ = TVr channels, which defines 
the dwell time through the system as tq/tq = M/2N. 
The dephasing time t^/tq = M/N3 is defined by the 
number A^3 of channels carried by the third, dephasing 
lead, and the Ehrcnfest time is given by 



A" 1 [lnA/ + 0(l)]. 



(85) 



2. The open kicked rotator coupled to an environment 

We extend the kicked rotator model to account for the 
coupling to external degrees of freedom. The exponential 
increase of memory size with number of particles forces us 
to focus on an environment modeled by a single chaotic 
particle. Therefore, we follow Ref. [68| and consider two 
coupled kicked rotators (with i = sys, env) 



J£ — H sys + H< 



H t = 



U 



(Pi - Pa) 



env 
2 



11. 



+ Ki cos(xi - x ) ^2 S(t - utq), 



sin(x 



sys -^env 



0.33)^<5(i-nro). (86) 



In this model, the interaction potential U acts at the 
same time as the kicks, which facilitates the construction 
of the S'-matrix. For this particular choice of interac- 
tion, the correlation length £ = L so that one expects 
a universal behavior of dephasing, as in Eq. (p}. The 
quantum representation of the coupled Hamiltonian is 
(82) a unitary (M sys M env ) x (M sys M env ) Floquet operator. 
We open the system (and not the environment) to two 
external leads by means of extended projectors Ptot = 
P' L ' (g) I onv + p( R ) g) I onv . A straightforward generaliza- 
tion of Eq. dHSJ) defines a (A^ T M cnv ) x (N T M cnv ) extended 
scattering matrix, from which we evaluate the conduc- 
tance via Eq. (|A11[) . We focus on the symmetric situation 
where the two leads carry the same number N of channels 
and average our data over a set of pure but random ini- 
tial environment density matrices ?7env(q, q';i = 0). We 
estimate the dephasing time from the e-induced broad- 
ening of two-particle levels in the corresponding closed 
two-particle kicked rotator, t^ 1 = 0.43(e/S e ff) 2 , with 
h c g = 2ir/M, the effective Planck constant. 



B. Weak- localization with dephasing 

To investigate weak-localization, we follow the same 
procedure as in Ref. [83] of taking a constant nonzero po, 



Ag 



Ag 




m n /m 

c 



Figure 8: a) Magnetoconductance curves Ag(xo) = g(xo) — 
g(0) for the double kicked rotator model (see text) with 
K syB = K cnv = 34.08 (Aeff « 2), r D /r = 8, £/L = 1 
and Hilbert space sizes M sys = 256, M em — 16. Different 
symbols correspond to different dephasing times t^/td = oo 
(H~gS — 0, circles), t^/td = 4.8 (h~gE ~ 0.25, squares), 
t^/td = 1.2 (fi~ff£ — 0.5, diamonds), t^/td = 0.3 (fr~g£ ^ 1, 
downward triangles) and t^/td = 0.07 (ft^e — 2, upward 
triangles) . Data are averaged over 25 different lead positions, 
each with 25 different quasi-energies and 10 different initial 
environment states, b) Magnetoconductance curves for the 
open kicked rotator with transparent dephasing lead (see text) 
with K = 34.08, td /to = 8 and Hilbert space size M = 256. 
Different symbols correspond to different dephasing times 
t>/td = oo (circles), t^/to = 5 (squares), t^/to = 1.25 (dia- 
monds), t^/td = 0.5 (downward triangles) and t^/td = 0.25 
(upward triangles). Data are averaged over 225 different lead 
positions, each with 50 different quasi-energies. 



while varying xq (which hence plays the role of a magnetic 
field). The obtained magnetoconductance in the absence 
of dephasing is Lorentzia n 45 i 83 , 



5 wl (0) 



1 + (m /m c y 



(87) 



with m c = 4n/K ^JMtd . In Fig. [5] we compare the sup- 
pression of weak-localization for the system-environment 
kicked rotator [panel a), top] and the dephasing lead 
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Figure 9: (Color online) Amplitude Sg(mo/m c — 2.4) — 
Sg(mo/m c — 0) of the weak localization correction to the con- 
ductance as a function of td /t^ for the double kicked rotator 
model (circles) and the open kicked rotator with transparent 
dephasing lead (squares). The black line gives the universal 
algebraic behavior (1+td /t { /,)~ 1 /4, and the red line is a guide 
to the eye, including an exponential decay with r^ 1 = 2.78, on 
top of the universal decay. 



kicked rotator [panel b), bottom]. For both models, we 
show five magnetoconductance curves, corresponding to 
five different ratios t^/td- All curves exhibit the ex- 
pected Lorentzian behavior vs. mo/m c , however, the am- 
plitude of the magnetoconductance is reduced as t^/td 
is reduced. This allows one to extract the r^-dependence 
of g wl . For the system-environment model, we found 
no significant departure from our analytical prediction, 
Eq. (pH) with = (since the interaction in Eq. (|8tj|) 
has £ = L). The same behavior is observed for the de- 
phasing lead model, as long as t^/td is large, however, 
compared to the system-environment kicked rotator, the 
damping of magnetoconductance accelerates as /te be- 
comes comparable to or smaller than one. When this 
regime is reached, magnetoconductance curves for the 
dephasing lead model lies significantly below those of the 
system-environment model, even when the latter corre- 
sponds to shorter dephasing times (compare in particular 
the upward triangles in both panels of Fig. [H]). This be- 
havior is further illustrated in Fig. [51 where we plot the 
amplitude 5g{m^/m c = 2.4) — Sg(mo/m c ) of the weak lo- 
calization corrections to the conductance as a function of 
td/t0. The data for the double kicked rotator nicely line 
up on the universal algebraic behavior (1 + td/t0)~ 1 /4 
without any fitting parameter. This is clearly not the 
case for the dephasing lead model, where an additional, 
exponential dependence on emerges. We attribute 
this to the exponential damping factor cx exp[— te/t^] 
of Eq. (f6"o) . The data presented in Figs. [5] and [H] for 
t^ 1 = exhibit a weak dependence on exp[— t^/to] only, 
with t£ < 0.8, which we attribute to terms of order one 
in Eq. ([8"5]). 

All collected data (including some that we do not 



present here) thus confirm qualitatively - if not quanti- 
tatively - the validity of Eq. ([55)) for the dephasing lead 
model. 



C. Conductance fluctuations in the dephasing lead 
model 

Conductance fluctuations were studied numerically in 
Ref. [23| for the dephasing lead model with a tunnel bar- 
rier of low transparency, and an exponential damping 
var(g) cx exp[— 2te/t^] was reported. Instead, here we 
consider a model in which the dephasing lead is trans- 
parently coupled to the system. In Fig. [TUl we show 
data for var(g) against t^/td, which is varied only by 
varying the width of the third, dephasing lead. There 
are four data sets (empty symbols) corresponding to a 
fixed classical configuration at different stages in the 
quantum-classical crossover, i.e. with increasing ratio 
M = L/X F G [128,8192]. As M increases, so does 
Te/t^,, however no change of behavior of var(g) is ob- 
served. These data are compared to a fifth set obtained 
in the universal regime, te/td *C 1, and the universal 
prediction of Eq. (|70|) (dashed line). These data confirm 
our analytical result, Eq. (|70|) . that var(g) exhibits no 
Ehrcnfcst time dependence for the dephasing lead model 
with perfectly transparent contacts. 



VI. CONCLUSIONS 

We have investigated the dephasing properties of open 
quantum chaotic system, focusing on the deep semiclas- 
sical limit where the Ehrenfest time is comparable to 
or larger than the dwell time through the system. We 
treated three models of dephasing. In the first one, the 
transport system is capacitively coupled to an external 
quantum chaotic system. For that model, we developed 
a new scattering formalism, based on an extended scat- 
tering matrix §, including the degrees of freedom of the 
environment. Transport properties are extracted from §, 
once the environment has been traced out properly. In 
that model, we find that, in addition to the universal al- 
gebraic suppression g wl oc (1 + td/t0) _1 with the dwell 
time td through the cavity and the dephasing rate t7 , 
weak-localization is exponentially suppressed by a factor 
oc exp[— tj/t^], with a new time scale depending on 
the correlation length of the coupling potential between 
the system and the environment. 

The second model we treated is that of dephasing due 
to a classical noise field. We show that the new time scale 
T£ plays the same role here as in the system-environment 
model. We then consider a classical Johnson-Nyquist 
noise model of electron-electron interactions. We show 
that £ ~ Af (and so t£ equals the Ehrcnfcst time) when 
dephasing is dominated by electron-electron interactions 
within the system, but that £ ~ D when dephasing is 




D 9 

Figure 10: Variance of the conductance vs. td/t^ for the open 
kicked rotator with K = 14. and td/to = 5 (empty symbols), 
transparently coupled to a dephasing lead. Different symbols 
correspond to different Hilbert space sizes (and hence different 
t£) M = 128 (squares, t^/td = 0.6), M = 512 (diamonds, 
te/td = 0.75), M = 2048 (upward triangles, t^/td = 0.9) 
and M = 8192 (downward triangles, t^'/to = 1.1). Addi- 
tional data for K = 144., r D = 25 and M = 2048 are also 
shown (full circles, r^ 1 /td = 0.08). The dashed line shows 
the universal behavior of Eq. ([2]). Unlike for weak- localization 
(see Fig. [8} and for the dephasing lead model with partial 
transparency^, the behavior of Sg 2 remains universal and 
shows no noticeable dependence on t^/to- Data are averaged 
over 50 different quasi-energies and from 50 (for N = 8192) 
to 500 (for N = 128 and 512) different lead positions. 



dominated by interactions between electrons in the sys- 
tem and those in a gate, a distance D away. 

The third model we treated is the dephasing lead 
model. We found a similar exponential suppression of 
weak-localization. To our surprise, however, it is the 
Fermi wavelength, not the dephasing-lcad's width, which 
plays a role similar to £ in that model. This inequivalence 
between the dephasing lead model and dephasing due to 
a real environment or classical noise can be most clearly 
seen in a situation where the interaction with a real en- 
vironment has £ ~ L. Then the environment induces 
only power-law dephasing, which is impossible to mimic 
with a dephasing lead. For the dephasing lead model, 
we also showed analytically and numerically that con- 
ductance fluctuations exhibit only power-law dephasing 
if the connection between the cavity and the dephasing 
lead is perfectly transparent. This is to be contrasted 
with the exponential dephasing observed for the dephas- 
ing lead model with tunnel barrier, and reflects the fact 
that the presence of tunnel barriers violates a sum rule 
otherwise preserving the universality of conductance fluc- 
tuations vs. te/td^. 

Related results have been obtained for conductance 
fluctuations in Ref. [24| . where different behaviors have 
been predicted for external sources - which give simi- 
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lar dephasing as our environment model - and internal 
source of dephasing - which qualitatively reproduce the 
prediction of Ref . [14[ . 

ACKNOWLEDGMENTS 

We thank I. Aleiner, P. Brouwer, M. Biittiker and M. 
Polianski for useful and stimulating discussions. While 



working on this project, C. Pctitjean was supported 
by the Swiss National Science Foundation, which also 
funded the visits to Geneva of P. Jacquod and R Whitney. 
P. Jacquod expresses his gratitude to M. Biittiker and the 
Department of Theoretical Physics at the University of 
Geneva for their hospitality. This work was initiated in 
the summer of 2006 at the Aspen Center for Physics. 



Appendix A: SCATTERING APPROACH TO TRANSPORT IN THE PRESENCE OF AN 

ENVIRONMENT 



Here we extend the scattering approach to transport to account for those environmental degrees of freedom which 
couple to the system being studied. We follow the lines of the derivation of the expression for noise presented in 
Ref. [8^ | . focusing on the two-terminal configuration. The following derivation is valid in the limit of pure dephasing, 
when there is no energy/momentum exchange between system and environment. 



1. Current operator and conductance 

In the presence of an environment, the current operator at time ( on a cross-section deep inside of lead a = L, R 
(where there is no system-environment interaction) reads 

L(t) = ~ jdE dE'e^ E - E '^ h ]T [at n (E')a an (E) - bi n {E')b an {E)] * W (Al) 



The second quantized operators and create and destroy incoming and outgoing system particles respectively. 
Since the environment particles carry no current, the current operator acts as the identity operator I env in the 
environment sub-space. We could write I en v in terms of second quantized operators, however for our present purpose, 
it is more convenient to write it as 

r dq|q)(q|. (A2) 

As in Ref. [84| we now back evolve the outgoing-states into incoming-states, this time with an 5*-matrix which also 
depends on the coordinates of the environment, 

b an (E){q\ = J2 / d qo S Q /3 TO (q,qo;£,£)a/y(£) (qol- (A3) 

Here, S Q( 3 ;n j (q, qo ; E, E) gives the transmission amplitude from channel j in lead f3 with energy E to channel n in 
lead a with energy E, while simultaneously, the environment evolves from q to q. We can set E = E, because 
throughout this article we only consider the regime of pure dephasing. Using (|A3[) we rewrite the current operator as 



/ a ( t ) = lJ dq J dEdE'e^ E - E '^ h £ 



(E) (*onv|q)(q|*onv) 



1,0 jk 

We next rewrite the first line of Eq. (|A4|) as 



dqo dq ^]^ (§a 7 ;ni(q,qi)) S aft „ fe (q,qo) a^E') ap k (E) (* cnv |q^) (q |* c 



(A4) 



dq S «« ( E ')^n (E) (* env I q) (q| * cnv ) 
n 

= j dqodq[ ) ^^ ( 5(q[ ) -qo)5 7Q 5 / 3 Q ^ fc a^( J B')a/3fe(-B)(*cnv|qo)<qo|*cnv). (A5) 



7,/3 jk 
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Finally, wc write the current operator in terms of the initial environment density-matrix ?/ cnv (qo, qd) = 

<qol*cnv)(*cnv|qo), 

I a (t) = J dq[ ) dq /( t rod )(q; ) ,q ;t)7 ?cnv (q ,q^), (A6) 
where we defined the reduced current operator as 

^ rcd) (qUo;i) = I /dq f dE&E'e^-z'W ^E^^^^^o^o)^^') a pk {E), (A7a) 

B j k (u,E,E';q' o ,q o ) = <5(q - q. )5 7a 6p a 5j k - /"dq^ (§a 7W (q, q )) §a/3 ; nfe(q, qo)- (A7b) 

The current is obtained by taking the expectation value of the current operator over the system, using 

((a\ .(£') a m (E))) = S J0 S jk S(E-E')f (E) (A8) 
where fp is the Fermi function in lead (3. The unitarity of § implies, 

J dq^ (s« 7i nj(qiCio)) §<5/3;nfc(q,qo) = <5(q„ -qo)^ 7 ^*- ( A9 ) 

S.n 

We use this latter equality to rewrite the Kronecker J's in Eq. (| A8|) . Finally the current in the left lead is 

« Ji » = ^ E / ^o^odq / d^(s £R; „ fe (q, qf,)) * § Lfl , ; „ fe (q, q ) [/l(£) - f R (E)] Vcav (q°, <&)■ ( Al 0) 

n.k 

In the limit of zero temperature in the leads, and assuming that the scattering matrix is not too strongly energy 
dependent, Eq. (|AlOj) leads to the linear conductance 

G = -r'Yl I dqodqodq(§i /fl; „ fe (q,qo)) Slr-m (q, qo)r?env(qo, qo), (All) 

with scattering matrices to be evaluated at the Fermi energy. From Eqs (jAlOj) and (jAllj) . we see that both current 
and conductance arc obtained by tracing over the environmental degrees of freedom of the square of the extended 
scattering matrix. Besides this prescription, these two equations arc extremely similar to their counterpart in the 
standard scattering approach to transport. We also note that conductance fluctuations can be obtained by squaring 
Eqs. ([A10]) and (|AlTj) . 

It is legitimate to expect that a complex environment - such as the chaotic system considered in this article - 
has a complicated initial wavefunction, which, under ensemble averaging, is uncorrelated with itself on all scales 
greater than the environment wavelength. This justifies us treating the initial environment state as (?7env(qo, qo)) ~ 
(2^) Jvd <5(q -q )/S onv . 



2. Current noise 



We follow similar steps as in the previous Section to calculate the zero-frequency current noise. However now we 
have two current operators, and hence two creation and two annihilation operators for the system^!. The environment 
has two I env operators, each of which we write in the form J dq|q)(q| to get the current time-correlator as 



((i (t)i a (o))) = 




(A12) 



The reduced current operator 1^ (q',q;i) is given in Eq. (|A7a[) . The zero-frequency noise power is obtained from 
the product of the deviations from the average current at times and t. Consequently it is proportional to 

J At J dqosdqordqor^/^qo^qoa;*) - <(i (rcd) (q i, q 03 ; *)» ) ( qoi; 0) - «/ (red) (q 3, qoi; 0))) ))) 

X»7cnv(q01,q l) ( A13 ) 
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There is only one trace over the environment here because we assume we measure the current as a function of time 
in a given experiment (with a given initial ?7 cnv ), and then extract the average current (and the deviations from it) 
from that data set. 

We need to take the following expectation value of products of creation/annihilation operators over the system^ 

((&l m (E 2 ) a-^E^al^iE,) ay n ,(E 3 ))) - «4 m (£ 2 ) M^M^'m't^) a, ( , n ,(E 3 ))) = 

<W V <W <W S(E 2 - E 3 ) 6{E 1 - E 4 )f a (E 2 )[l T fp(Eh.)], (A14) 

where the minus (plus) sign stands for fcrmions (bosons). From this we finally get the zero-frequency noise power 
Saaf(0) = j J dE /dqo 3 dq 01 dq^ 1 ^^^7(a,i?,i?;q; )1 ,qo3)S^r(a',i?,i?;qo3,qoi) 

x(fp(E)[lTf-y(E)] + [lTfa(E)]fp(E)) ?7cnv(qoi,qoi). (A15) 



All the relevant information for shot noise is contained in the diagonal S aa , a = L,R. Shot noise is obtained by 
calculating this latter expression in the limit of low temperature but finite voltage bias V between the two leads. In 
that case, it is easily checked that the contribution to B arising from the first term on the right-hand side of Eq. (|A7bj) 
docs not contribute, and one gets 



r r t 

S aa (0) = / dE ^2 / dqsdqidqogdqoidq^^] ^ (&a/3;m'm(<ia, Qoi; E )j ^a^m'ni^, qo3; E) 

/3^7 m,nm',n' 

x(§ Q7; „'„(qi,qo3;^)) t §aft„' m (qi,qoi;^) (fp(E)[l T f 7 (E)} + [1 T f a (E)}fp(E^rj cnv (q 01 ,^ 01 ).(Al6) 

We finally assume a slow dependence of S on E, in which case the integral over the energy is easily performed, giving 
a factor eV. For a two- lead device (L,R), we find that 

Srr(O) = '—j^- J dqsdqidqosdqoidqo! ^ (§RL;m> m (q.3, Qoi 5 E fj ^ $ R R ;m >n(q.3, qo35 E) 

xfSi?fl. n / n (qi,qo3;-B)) Sfli ;n 'm(qi, qoi; E) r? en v(qoi, Qoi)- (A17) 
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